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Semiclassical periodic orbit theory is used in many branches of physics. However, most applications 
of the theory have been to systems which involve only single particle dynamics. In this work, we 
develop a semiclassical formalism to describe the density of states for two noninteracting particles. 
This includes accounting properly for the particle exchange symmetry. As concrete examples, we 
study two identical particles in a disk and in a cardioid. In each case, we demonstrate that the 
semiclassical formalism correctly reproduces the quantum densities. 



I. INTRODUCTION 

Semiclassical physics has experienced a resurgence of in- 
terest, largely due to the work of Gutzwiller [Q, Balian 
and Bloch ||] and Berry and Tabor ||^. (For recent re- 
views see pH-) These works showed that if we separate 
the density of states into smooth and oscillatory com- 
ponents, then the oscillatory part is related to the dy- 
namics of the underlying classical system via periodic or- 
bits. This complements the earlier work of Weyl, Wigner, 
Kirkwood and others who showed that the smooth com- 
ponent is related to the geometry of the classical phase 
space. Actually, the two components are related in a sub- 
tle way 1^,01 since the complete geometry imparts the full 
dynamics and vice-versa. 

Most of the theoretical work has concentrated on the 
single particle density of states, however, there are some 
notable exceptions, namely [H^^Oj. In ||^ the focus is on 
the average level density and its extension to systems of 
identical particles. Specifically, the authors consider a 
system of N fcrmions in one dimension. Their Weyl for- 
mula for fermions works well for attractive two body in- 
teractions, but overestimates the quantum staircase func- 
tion when there are repulsive two body interactions. The 
author of |^ develops a generalization of the canonical 
periodic orbit sum for the special case of N interacting 
spinless fermions in one dimension. It is assumed the 
periodic orbits are isolated and therefore it is most ap- 
plicable to fully chaotic systems. The author also con- 
siders a system of noninteracting fermions and writes the 
many body level density as a convolution integral involv- 
ing one body level densities. Finally, we mention ||ic| ] 
which presents an expansion of the periodic orbit sum in 
terms of the particle number using ideas from |^,^ . 

Similarly, most of the applications of semiclassical the- 
ory have been to systems which involve only single par- 
ticle dynamics. Here, we mention some exceptions. The 
authors of [0 extend the study of scars |12| to classically 
chaotic few body systems of identical particles. A study 
of the eigenfunctions of an interacting two particle sys- 
tem can be found in [ p^ . The semiclassical approach to 
the helium atom, which can be understood as two inter- 



acting electrons in the presence of a helium nucleus, has 
been studied in . We also mention the novel applica- 
tions of semiclassical theory to mesoscopic physics 
For example, orbital magnetism has been studied semi- 
classically for diffusive systems in llfl] and for ballistic 
systems in |l^Jl^. For reviews see |19[| . 

Ultimately, one would like to study an arbitrary num- 
ber of interacting particles in any kind of potential. In 
the present work, we begin by exploring the structure of 
the trace formula for two noninteracting particles includ- 
ing an examination of the decomposition into bosonic 
and fermionic spaces. This sets the stage for the inter- 
acting A'^— body problem to be explored in a subsequent 
publication |2^. The method employed here uses the 
fact that the two particle density of states is the auto- 
convolution of the single particle density of states. Sub- 
sequently, we decompose the semiclassical two particle 
density of states into three distinct contributions, each 
of which corresponds to specific dynamical properties of 
the system. Of particular interest is the contribution 
which corresponds to two particle dynamics. 

Billiards have served as prominent model systems in 
quantum chaos. They combine conceptual simplicity (the 
model of a free particle in a box) while allowing the full 
range of classical dynamics, from integrable to chaotic. 
Therefore, as initial applications of the formalism, we 
study two noninteracting identical particles in a disk and 
in a cardioid. The former problem is integrable while 
the second is chaotic so these two examples provide a 
direct test of the formalism in the two limiting cases of 
classical motion. In both cases, we find the semiclassical 
formalism does a good job of reproducing the quantum 
density of states. 



II. BACKGROUND THEORY 

A. Single Particle Semiclassical Theory 

In this section, we review the formalism for the semi- 
classical decomposition of the single particle density of 
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states. Let {ei} be the single-particle energies so that 
the single particle density of states is 

Pi(6) = ^<5(e-e0, (1) 

i 

where the subscript 1 indicates that it is a single particle 
density. A fundamental property of the quantum density 
of states is that it can be exactly decomposed into an 
average smooth part and an oscillatory part ^ 

pi(e) = pi(e)+pi(e). (2) 

There are various approaches for calculating these quan- 
tities 1^. For example, in systems with analytic poten- 
tials, the smooth part may be obtained from an extended 
Thomas- Fermi calculation which is an asymptotic expan- 
sion in powers oi h. In billiard systems, where the par- 
ticle is confined to a spatial domain by the presence of 
infinitely steep potential walls, the smooth part may be 
obtained from the Weyl expansion. In two dimensional 
billiards with piecewise smooth boundaries and Dirichlet 
boundary conditions, the first three terms of the Weyl 
expansion is |^ 

+ (3) 

where a = 2m /h^, A is the area, £ is the perimeter and 
JC = J- I dlK{l) + —Y ^^—^ (4) 

is the average curvature integrated along the boundary 
with corrections due to corners with angles 6i. The oscil- 
lating part is obtained from semiclassical periodic orbit 
theory, and in particular the various trace formulas for 
pi{e) of the form [| 

^i(^) « ^E^r(e)cos Q5r(6) - ar|) . (5) 

F denotes topologically distinct periodic orbits, 5'r(e) is 
the classical action integral along the orbit F. The am- 
plitude Ar{e) depends on energy, the period of the cor- 
responding primitive orbit, the stability of the orbit, and 
whether it is isolated or non-isolated. The index err de- 
pends on the topological properties of each orbit. For 
isolated orbits, it is just the Maslov index. For noniso- 
lated orbits, there may be additional phase factors in the 
form of odd multiples of 7r/4 which we account for, in a 
slight abuse of notation, by allowing or to be half-integer. 
In the case of non-isolated orbits, F denotes distinct fam- 
ilies of degenerate orbits. The amplitude of an isolated 
orbit is given by the Gutzwiller trace formula jl] 

Ar{e) = (6) 
,/ det(Mr - /) 



where T^{e) is the period of the primitive orbit 7, corre- 
sponding to F (i.e. F is an integer repetition of 7) and 
Mr is the stability matrix of that orbit. 

B. Quantum Two Particle Density of States 

Now suppose we have a system of two identical nonin- 
teracting particles. The total Hamiltonian is the sum of 
the single particle Hamiltonians and it follows that the 
energies of the composite system are just the sums of the 
single particle energies. The analogue of (Q) is then 

P2(£;)=^<5(i?-(e, + 6,)). (7) 

hi 

A useful relation is that the two particle density of states 
is the autoconvolution of the single particle density of 
states: 

P2iE)^ [ depiie)piiE-e)^pi^piiE), (8) 
Jo 

as can be verified by direct substitution. In fact, this 
works even if the particles are not identical, where the 
full density is still the convolution of the two distinct 
single-particle densities. This would also apply to a single 
particle in a separable potential, which is mathematically 
equivalent. Rather than encumber the notation to explic- 
itly allow for this possibility, we defer this discussion to 
Appendix A, where some formulas for nonidentical, non- 
interacting particles are presented. 

We can decompose the two particle density of states 
for a system of two identical particles into a symmetric 
and an antisymmetric density, 

P2iE)=psiE)+pAiE). (9) 

We shall use the terms symmetric/antisymmetric and 
bosonic/fermionic interchangeably. Each partial density 
may be obtained using a projection operator onto the 
relevant subspaces resulting in 

PsM(i^) = ^P2(£;)±ipi(|)). (10) 

We seek semiclassical approximations to these quantum 
expressions, a topic which is pursued in the following 
sections. 

III. SEMICLASSICAL CALCULATIONS FOR 
THE TWO PARTICLE SYSTEM 

Decomposing the single particle density into its smooth 
and oscillatory components as in (^) gives a decomposi- 
tion of the two particle density of states into three dis- 
tinct contributions. 
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pI%E) = pi * piiE) + 2pi * pi{E) + pi * pi{E). (11) 

The first term is a smooth function of energy since the 
convolution of two smooth functions results in a smooth 
function. This is followed by a cross term and finally by 
a purely oscillating term. The cross term is also an os- 
cillating function. At first, this may seem incorrect since 
the convolution of a smooth function with an oscillating 
function usually yields a smooth function. As we will 
show, the oscillatory nature of the cross term is due to 
contributions from the end-points of the convolution in- 
tegral. Physically, the smooth term does not depend on 
dynamics since it corresponds to the Weyl formula in the 
full two-particle space. The cross term depends only on 
single particle dynamics because it corresponds to the 
situation where one particle is stationary and the other 
particle is evolving dynamically on a periodic orbit. It is 
only the last term which contains two particle dynamics 
in the sense that both particles are evolving dynamically 
on periodic orbits. Hence, we will refer to the last term 
as the dynamical term. 

We find a general expression for pi * pi{E) by sub- 
stituting a generalized trace formula for pi (E) and then 
evaluating the resulting convolution integral using the 
method of stationary phase. Using (^, the dynamical 
term can be written as 



Pi * Pi{ 



^ ' ri,r2 ° 



cos Q^r, (e) - ar, |) cos Q^r, {E ^ e) ^ cjr, |) . (12) 

To evaluate this asymptotically, we should include all 
critical points in the integration domain. Specifically, 
this integral has a stationary phase point within the in- 
tegration domain and finite valued endpoints. We shall 
show that the stationary phase point corresponds to the 
situation where both particles are evolving dynamically 
with the energy partitioned between the two particles in 
a prescribed way. The endpoint contributions must be 
evaluated at energies such that one of the particles has 
all of the energy while the other has no energy. However, 
this contradicts our assumption that both particles are 
evolving — this is the definition of the dynamical term. 
Moreover, if we were to evaluate this contribution, the 
result would be meaningless since it involves using the 
trace formula at zero energy where it is known to fail. So 
we shall omit the contributions from the endpoints; this 



is discussed more fully in section IV D and in Appendix 
B, as well as in reference p2| . 

Hence, we evaluate the integral in (^2|) using only the 
stationary phase point. To leading order, we can extend 
the integration limits over an infinite domain. Writing 
the cosine functions as complex exponentials yields four 
integrals; the first is 



deAr, (e)Ar, {E ~ e) exp ( - (St, (e) + St, (E - e)) 



AtAEo)At,{E- Eo). 



2-Kh 



|T(ri,r2,£;)| 



e^'p\'-{STAEo) + ST,{E-Ea))+iv^ 



(13) 



where 



T(ri,r2,i?) = 



O^StA^) , d^ST^E-e] 



Eo 



z/-sign(T(ri,r2,S)). 
Eq is determined from the stationary phase condition 
'5S'ri(e) dST,{E-£ 



(14) 



de 



de 







Tt,{E^) = Tt,{E-Eo) (15) 



where we have used the fact that the derivative of the 
action with respect to energy is the period. Eq is the 
energy of particle 1, E — Eq is the energy of particle 2 and 
E is the total energy of the composite system. The saddle 
energy Eq has a precise physical interpretation; Eq.(|l5|) 
says that the energies of the two particles are partitioned 
so that the periods of both periodic orbits are the same. 
In other words, at Eq, we have orbits which are periodic 
in the full two particle phase space since after the period 
T both particles return to their initial conditions. 

The next integral has the same stationary phase con- 
dition as the first integral and is its complex conjugate. 
The third integral is 



deApi (e) (E - e) exp 



(^r.(e) 



SrAE~e)) 
(16) 



and has no stationary phase point since setting the first 
derivative of the action to zero yields the stationary phase 
condition 



TtAEo) = -Tt,{E-Eo). 



(17) 



The trace formula only involves orbits with positive pe- 
riod, so we ignore this possibility. The last integral is the 
complex conjugate of the third and will also be ignored. 

Adding the contributions from the first two integrals, 
we arrive at the two particle trace formula: 



pi * pi{E) 



(27r;i)3/2 



E 

ri.ra 



AT,{En)AT,{E - Eq) 

V|T(ri,r2,i?)| 



1 TT 

COS I -{StAEo) + St^E ~ Eo)) - {(JT, +CTrJ2 



This result possesses the intuitive properties that, other 
than factors arising from the stationary phase analysis. 



(18) 
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the actions and Maslov indices are additive and the am- 
phtudes are multiphcative. We note that this saddle- 
point analysis fails for the simplest problem in physics, 
the harmonic oscillator, where T = 0. This is because 
the two-particle harmonic oscillator has a higher degree 
of symmetry than we are accounting for here. This is a 
nongeneric property specific to the harmonic oscillator. 
We also stress that we have made no assumption about 
the stability or structure of the orbits. They can be iso- 
lated, stable or unstable or come in families. There are 
also problems with coexisting isolated orbits and families, 
such as those of the equilateral triangle billiard |^,|| . 
Note that the overall h dependence is not multiplica- 

1 /2 

tive but picks up an additional factor of h ' from the 
stationary phase integral. For isolated orbits, the ampli- 
tudes A are independent of Ti and the expression (|lj) has 
a prefactor as opposed to the 1/h in the amplitude 

of the single particle trace formula. The fact that the Ti 
dependence is different implies that the periodic orbits 
of the full system come in continuous degenerate families 
rather than isolated trajectories, which in turn implies 
that there exists a continuous symmetry in the problem 
Ip^ . This is an important point which we will address 
in a companion paper |2^. (It was also noted in 0.) 
Nonetheless, it may be helpful to give a brief explana- 
tion here. Imagine the full phase space periodic orbit F 
consists of particle 1 on a periodic orbit Fi with energy 
Eq and particle 2 on a distinct periodic orbit F2 with 
energy E — Eq. We can define t = to be when particle 
2 is at some prescribed point on F2. Keeping particle 2 
fixed, we can change the position of particle 1 on Fi to 
generate the initial condition of a distinct but congru- 
ent periodic orbit in the full phase space. Continuous 
time translation of the initial condition on Fi generates 
a continuous family of congruent periodic orbits in the 
full phase space. Since the time translational symmetry 
can be characterized by a single independent symmetry 
parameter, the Ti dependence is O stronger than 

for a system with isolated periodic orbits p^,p[|. 



IV. TWO PARTICLE QUANTUM BILLIARDS 

As an application of the formalism developed in section 
III, we consider the quantum billiard problem. Billiards 
are two dimensional enclosures that constrain the mo- 
tion of a free particle. Classically, a particle has elastic 
collisions with the walls and depending on the geometric 
properties of the domain, the dynamics are either regu- 
lar or chaotic. For the nonintcracting problem, the two 
particles move independently of each other. In a billiard 
system, classical orbits possess simple scaling properties. 
For instance, the action of an orbit F, S'r(e) — \^2meLr 
and the period of the orbit is 



Trie) 



dS- 



re 



'2m Lt 



de 



(19) 



The parameter a = 2m /h^ already appeared in Eq. (|^); 
it will recur often. For example, in all final expressions, 
the energy occurs with a; this is a result of the scaling 
property (the quantity aE having the units of 1/length^.) 
In the theoretical development, it will be convenient to 
retain a and use it to keep track of relative orders in the 
semiclassical expansions (since it contains K). However, 
once we have the final expressions, we are free to set it 
to unity for the purposes of numerical comparisons. 

We also mention that for billiards, it is common to ex- 
press the density of states in terms of the wave number 
/c, where e = fc^/a so that p{k) — 2kp{e)/a. This is con- 
venient since k is conjugate to the periodic orbit lengths 
L. Therefore, many of our results will be quoted as a 
function of fc, although it should be stressed that all con- 
volution integrals must be done in the energy domain. 
Thus, we shall write 

pf{k) = (pi * pi) (fc) + 2 (pi * pi) (fc) + (pi * pi) (fc). 

(20) 

Here, it is understood that each of the functions in brack- 
ets is first evaluated in the energy domain and then con- 
verted to the k domain through the Jacobian relation 
above. This will always be the case when the argument 
is fc, so that we will not always write brackets around 
the various functions. In terms of the wavenumber k, the 
decomposition (noh becomes 



PS,A{k)^\[p,{k)±l=p, (^_| 



(21) 



A. Smooth Term 



The smooth part is defined by the convolution integral 



pi * pi{E) = i depi(e)pi(£;- 
Jo 



(22) 



where pi is given by the Weyl expansion. The expansion 
in (H) is taken only to order Ti". Hence, after expanding 
the integrand in (^) , it is formally meaningless to include 
terms that are O (1/fi.) since there are corrections of the 
same order in Ti that have not been calculated. Ignoring 
these terms and performing the necessary integrations, 
the smooth term is found to be 



pi * pi{E) 



167r2 



-E 



87r2 



E^ 



a£2 aAlC 



647r 



27r 



(23) 
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B. Cross Term 

We next convolve pi term by term with pi. Asymp- 
totically, each convolution integral receives contributions 
from the upper and lower endpoints. However, we shall 
only include one of these, namely the endpoint for which 
the trace formula is not evaluated at zero energy. As in 
section [I], we neglect the other endpoint for reasons ex- 
plained in IV D and Ref. |23]. This is also discussed in 



Appendix B, where we evaluate the various integrals for 
the cross term exactly using isolated billiard orbits and 
show explicitly that an appropriate asymptotic expansion 
of the exact expression leads to consistent results. 

After convolution, we find the area term involves the 
integral 



Rej^ deAr{E~e) exp (i^/a{E 



- e)Lr - *o-r-^^ 



(24) 



The lower endpoint e = corresponds to the physically 
meaningful situation while the upper endpoint is spuri- 
ous in the sense mentioned above and discussed in detail 
below. Hence, to leading order, we can remove the am- 
plitude factor from inside the integral, Taylor expand the 
argument of the exponential and extend the upper limit 
to infinity. This leads to 



aA Ar 



r ^ 



^ir-ar|-|). (25) 



By similar logic, the perimeter term and curvature terms 
are 




lie (E) ~ ^ cos (VaELr - err ^ 



Note that all amplitudes and periods in (g5|) and (|2q) 
are evaluated at the system energy E. Recall a oc 
so that after convolution the sequence is an expansion in 
powers of \fh and not in powers of h as for the original 
Weyl series (||). We also note that the first correction to 
Ij^ may be of the same order as (as happens for the 
disk [^) and should be included if this is the case. We 
then have p * p ^ Ia + Ic + lie- 



C. Dynamical Term 

In this section, we derive a general expression for the 
dynamical term that is valid for any billiard problem. To 
this end, the first task is to determine the saddle energy 
from the stationary phase condition. Inserting ( |T9| ) into 
M) yields 



which implies 

Ea 
E 



7^2 



/Eq \JE — Eq 



E — En 



7^2 



E 



and 



T(ri,r2,i?) = 



'2m. 



i 2 



,5/2 



4^;3/2 



i 1 i 2 



(27) 



(28) 



(29) 



Clearly = — 1. We then substitute these results into 
( |l8| ) to obtain the two particle trace formula for billiards 



pi * piiE) 



v^ai/4(27r?i)3/2 




ArAEo)ArAE - Eo) 



If the single particle periodic orbits are not isolated, then 
one must make direct use of the corresponding single par- 
ticle amplitudes in (^ evaluated at the appropriate en- 
ergies. We will show an explicit example of this when we 
analyze the disk billiard. Note the amplitudes Ar typi- 
cally have an energy dependence so one cannot make any 
general statements about the energy dependence of this 
term except that the greater the dimensionality of the 
periodic orbit families, the greater the energy prefactor. 
For example, for the disk, it turns out to be E^^'^. 

If the single particle periodic orbits are isolated, the 
amplitudes are given by (0), which for billiards is 



Ar(e) 



yfah 



det [Mr -I 



(31) 



In this case, the Gutzwiller amplitudes are evaluated at 
Eq and E — Eq, so we again make use of (|2^). After some 
algebra and simplification, we find 



Pi * Pi{E) 



a 



3/4 



(27r)3/2£;i/4 



E 

cos 



det (^Afn - /) dot (^Mr^ - /) 



Note the E''-^/^ prefactor which implies that the ampli- 
tude decays weakly with energy. This is the same prefac- 
tor that occurs in the single particle disk problem. This 
is not a coincidence, but arises from the fact that in both 



(30) 



(32) 
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problems the periodic orbits come in one parameter fam- 
ilies. Also, one must be careful to distinguish between 
Lr, the length of a periodic orbit and i-y, the length 
of the corresponding primitive periodic orbit. In general 
Lr = nrL^ where nr is the repetition index of that orbit. 

D. Spurious Endpoint Contributions 

As mentioned above, when confronted with convolution 
integrals, it is natural to analyse them asymptotically. 
This involves identifying the critical points and doing 
appropriate expansions in their neighbourhoods. In our 
work, these critical points are either stationary phase 
points or endpoints. The power of semiclassical meth- 
ods is that each critical point can be given an immedi- 
ate physical interpretation. For example, the stationary 
phase point in the dynamical term is found to be that 
energy such that the two particles have the same period 
so that the motion is periodic in the full two-particle 
phase space. This is intuitively reasonable. However, the 
same integral also has endpoints with finite valued con- 
tributions. We could do an asymptotic calculation in the 
vicinity of these points, but we can argue immediately 
that the result is spurious and not physically meaning- 
ful. 

Recall the trace formulas are asymptotic in h which 
typically also means asymptotic in energy. At the end- 
points, one of the trace formulas is evaluated at small 
energy where it is known to be invalid. Alternatively, we 
can substitute for the trace formula any expression which 
is asymptotically equivalent to it and expect all meaning- 
ful results to be invariant to leading order. If we do this, 
we will find the endpoint contribution changes while the 
stationary phase contribution remains invariant, to lead- 
ing order. 

A further argument is that the structure of the end- 
point contribution will be incorrect. Typically, it will be 
a sinusoid with an argument which does not depend on 
energy, but only depends on the properties of one of the 
orbits. Hence, it will have the same asymptotic structure 
as the cross term. However, we know that the cross term 
completely describes all such contributions and any fur- 
ther contribution with the same structure must be spu- 
rious. 

Similarly, when we evaluate the cross term, we have 
two endpoint contributions. At one of these, we are eval- 
uating the trace formula at some finite energy, which is 
reasonable. This endpoint corresponds to orbits which 
are periodic in the full phase space and in which one par- 
ticle evolves on a single particle periodic orbit with all 
the energy, while the other remains fixed at some point in 
phase space with zero energy. At the other endpoint, we 
are evaluating the trace formula at zero energy, which is 
problematic. This corresponds to the contradictory situ- 
ation in which the evolving particle has zero energy while 



the fixed particle has all the energy. In addition, upon 
inspection of this endpoint contribution, we find a func- 
tion which is not oscillatory in energy and therefore has 
the same asymptotic structure as the smooth term. How- 
ever, the smooth term already completely describes the 
average behaviour of the two particle density of states 
and any further contributions with the same structure 
must be spurious. 

These situations are further examples of a general sit- 
uation described in Ref. |Q where it was shown that 
when integrating over the trace formula to obtain phys- 
ical quantities, one should include all critical points ex- 
cept ones at which the trace formula is evaluated at zero 
energy. Such contributions should simply be ignored as 
spurious. In |22| , the application was to thermodynamic 
calculations, but the principle is precisely the same. In 
Appendix B, we show the result of evaluating the cross 
term exactly for isolated orbits. An asymptotic analysis 
of this result leads to two terms which we can identify as 
coming from the two endpoints. One has the form used 
in this paper while the other is clearly spurious. 

V. TWO PARTICLE DISK BILLIARD 

In this section, we apply our results to the problem of 
two identical noninteracting particles moving in a two 
dimensional disk billiard of radius R. Quantum mechan- 
ically, this problem is a simple extension of the one body 
problem. Nevertheless, the spectrum has some interest- 
ing features which we discuss below. 

A. Quantum Mechanics 

For the disk billiard, a general two particle state can be 
written as 

\mi ni,TO2 ^2) = |mi ni) (g) |to2 712) (33) 

where the azimuthal quantum numbers mi, TO2 — 

0, ±1, ±2, . . . and the radial quantum numbers ni, n2 — 

1, 2, 3, . . .. We shall also use a more compact notation 

N2) = \mi Til, 1712 ^2) where denotes a pair of in- 
tegers (m,n). We can immediately write down the wave 
numbers of the two particle system as 




where Zn denotes the nth zero of the mth Bessel func- 
tion Jm{z). The set of all two particle states is given by 
{\Ni,N2)}. 

The spectrum is highly degenerate. A typical state 
|A^i, A^2) is 8- fold degenerate since we can reverse the sign 
of either mi or m2 or interchange the two particles and 
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the resultant state has the same energy. However, if ei- 
ther mi or TO2 is zero or if A^i = 7V2 then the state is 
4-fold degenerate. If toi = m2 = and Ni ^ N2, then 
the state is 2-fold degenerate whereas if mi = m2 = and 
Ni = N2, then the state is nondegenerate. If the particles 
are in distinct states, the degenerate multiplets divide 
evenly between the symmetric and antisymmetric spaces. 
However, if the particles are in the same state, A^i = N2, 
it is somewhat less trivial. If A'^i = N2 and mi ~ m2 7^ 0, 
there is a 4- fold degenerate set of states: \m n,m n), 
\—m n, — m n), |m n, — m n) and \—m n,m n). The first 
two states belong to the symmetric space. From the sec- 
ond two states, we can construct one symmetric and one 
antisymmetric combination. (This is analogous to cou- 
pling two spin 1/2 states to construct a 3-fold symmetric 
S — \ state and a nondegenerate antisymmetric S" = 
state.) If Ni = N2 and mi = m2 = 0, this yields the 
state |0n. On), which is singly degenerate and belongs 
to the symmetric space. 



B. Semiclassical Density of States 

We first review the semiclassical decomposition of the 
single particle density of states. The smooth part of the 
density of states may be obtained using the general result 
for two dimensional billiards (H). In fact, many higher 
order terms have been calculated 0. But, for our pur- 
poses, it suffices to use the first three terms as in (^) with 
A = 7ri?2, C = 2ttR and K. = 1/6. 

The oscillating part of the level density can be obtained 
using trace formulas for systems with degenerate families 
of orbits. The periodic orbit families may be uniquely 
labelled by two integers (v, w) where v is the number of 
vertices and w is the winding number around the center. 
The two integers must satisfy the relation v > 2w. The 
length of an orbit with vertex number v and winding 
number w is given by Lvw = 2vi?sin (ttw/v). With this 
notation, the trace formula for the oscillating part of the 
density of states is [E6[ 
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FIG. 1. (Top) The quantum density of states for two iden- 
tical particles in the disk billiard. (Middle) Bosonic density 
of states. (Bottom) Fermionic density of states. In each case, 
the heights indicate the degeneracy of the state. 



The quantum density of states 

P2{k)= ^ 5{k-kMiN2) 

Ni,N2 



(35) 



and the corresponding symmetric and antisymmetric 
densities are shown in Fig. ^ as a function of the 
wavenumber k. Note that in this figure some of the peaks 
have different degeneracies in the symmetric and anti- 
symmetric densities, as discussed above. 



v3/4 



3/2 



2\/2^ei/4 ^ v2 




vr TT 

^^2 + 4 



(36) 



where the sum goes from w = 1 • • • 00 and v = 2w • • • 00 
and the degeneracy factor I?vw, which accounts for neg- 
ative windings, is 1 for v = 2w and 2 for v > 2w. Com- 
paring (Bra) with the general form (ph, we identify 



^vw(£) 



27ra3/4^P r3/2 



3v 



4v2ei/4 
1 

2 



(37) 



Adding the smooth and oscillating terms gives the semi- 
classical approximation to the single particle density of 
states which we denote by Pi'^(e). 

To evaluate the semiclassical approximation to the two 
particle density of states, we must evaluate the smooth, 
cross and dynamical terms. The smooth term can be 
taken from Eq. (p3[) to be 



3/2^3 



E 



37r + 4 



(38) 



The arguments of the previous section and in particular 
Eqs. (psi ) and (|2^ ) lead to the cross term 

Pl * Pl{E) « ;-^= 2^ ■ 



4V27r 



-y^xvw cos (-i>™ - ^) 



VW 

7r\ /I 



Xvw COS -S^ 



(39) 



7 



where <J)vw = VaEL^^ - 3v7r/2 + 7r/4 and Xvw = 
\/L^ / {aEy/ ^ R. We have also included the first cor- 
rection to the area term, Ia{E) which appears in the 
third term above [ p5[ |. The dynamical term can be ob- 
tained using (|30|). Noting that Ti in ( |30| ) corresponds to 
the pair of integers (vi,Wj), the result is 



pi * pi{E) 



,5/4^1/4 



-3/2 



4\/27r ^ V-'-^ v: 

* ViWi,V2W2 \i = l 



cos (y^Li2 - 3 (vi + V2) I + ^) (40) 
where we defined Li = L^ ^i, T)i = Wi and L12 = 




FIG. 2. The cross term (|39| ) of the semiclassical density of 
states for two identical particles in the disk billiard. In this 



case, we truncate the sum in 



at Wmax = 50, Vmax = 100. 



The circles indicate the level sequence of the one body prob- 
lem obtained from EBK quantization. Note the kinks that 
occur at these positions. 



C. Numerics 

For numerical purposes, we take a = h = 1 and i? = 1 so 
the single-particle energies are just the squares of the ze- 
ros of Bessel functions. Since we can only include a finite 
number of orbits, the periodic orbit sums must be trun- 
cated. As a representative case, we truncate the sum 
in (|3^) at w^ax — 50,Vniax = 100 (see Fig. 0) and use 



the same limits to truncate the quadruple sum in (40) 



This is a relatively small set of orbits, yet it does very 
well in reproducing the peaks of the quantum density of 
states. As an illustration, we show the first few peaks of 
( pO| ) in Fig. 1^. We calculated the semiclassical density of 
states (po[) on the interval < fc < 11. After doing so. 



we found only two sets of two peaks which were not re- 
solved. These are shown in Fig. ^. Obviously, using more 
orbits will produce better results, but this increases the 
computation time because of the quadruple sum in (|4^). 
(Although one can reduce the computational overhead 
by limiting the sum to orbits whose amplitude exceeds 
some prescribed threshold [p5|). 
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FIG. 3. (Top) The first few peaks of the semiclassical den- 
sity of states (EOl). (Middle) Semiclassical approximation to 
the bosonic density of states. (Bottom) Semiclassical approx- 
imation to the fermionic density of states. In each case, the 
circles and stars represent the appropriate level sequences ob- 
tained from EBK quantization and quantum mechanics, re- 
spectively. Note the positions of the peaks more closely re- 
produce the EBK spectrum. 

As an additional test, we want to determine whether 
( pO| ) gives the correct degeneracies. We could do this 
by integrating the area under each of the peaks. How- 
ever, a simpler procedure is to do a Gaussian smoothing 
by convolving p^{k) with an unnormalized Gaussian of 
variance a: 



pf{k)*G, 



= r 

Jo 



Ak'pf{k')G^{k-k') (41) 



where 



G„{k) = exp(-fcV2cr2). 



(42) 



and a is the smoothing width. The reason for this is 
that if the variance a of the Gaussian is larger than the 
intrinsic width of a peak in the semiclassical spectrum, 
then each peak acts like d5{k — fc„) with respect to the 
Gaussian. Thus, the integral in ( pl| ) becomes dGa{k — 
kn) OT d a,t k = kn- Of course, this is invalid when the 
spacing between two adjacent peaks is smaller than about 
a. Some examples are discussed in the next section. 
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FIG. 4. Two sets of two peaks in the semiclassical density 
of states ( po| ) that are not resolved. (Left) The middle peak 
is not resolved into the two peaks at fc = 7.4163 (correspond- 
ing to the quartet {|0 1,±1 2) , |±1 2,0 1)}) and k = 7.4423 
(corresponding to the octet {|±1 1, ±3 1) , |±3 1, ±1 1)}). 
The corresponding EBK quartet and octet energies occur 
at = 7.3831 and k = 7.3932 respectively. Note that 
AksBK ~ 0.0101 and AkqM = 0.026 so that the spac- 
ing of the two unresolved levels is smaller in the semi- 
classical spectrum than in the quantum mechanical spec- 
trum. (Right) The middle peak is not resolved into the 
two peaks at = 9.4641 (corresponding to the quartet 
{|±1 1,0 3) , |0 3, ±1 1)}) and k = 9.4829 (corresponding 
to the octet {|±3 1, ±1 2) , |±1 2, ±3 1)}.) The correspond- 
ing EBK quartet and octet energies occur at = 9.4359 
and k = 9.4456 respectively. Here, AksBK = 0.0097 and 
AkQM = 0.0188. 

We also studied the symmetrised densities by using the 
expression (21) for both the quantum and semiclassical 
densities and convolving as above. The periodic orbit 
sums in the oscillating parts of the one and two body 
densities were truncated in the standard manner as be- 
fore. The result of this numerical procedure is shown in 
Fig. ^. Clearly, the semiclassical approximations repro- 
duce the correct degeneracies of the quantum spectrum 
as well as the approximate positions. 



FIG. 5. (Top) The smoothed semiclassical densitvof states 
obtained from numerical convolution of ( po| ) with (ji^). Note 
the artifacts of the single particle EBK spectrum which occur 
at the positions marked by an "X". (Middle) The smoothed 
semiclassical bosonic density of states obtained from numeri- 
cal convolution of (^l|) (with the -I- sign) and (^^. (Bottom) 
The smoothed semiclassical fermionic density of states ob- 
tained from numerical convolution of ( pl| ) (with the - sign) 
and (^^. In each case the sequence of dots represent the 
corresponding EBK spectrum and a — 0.0125. 



D. Discussion 

In 1^ , it was noted that the trace formula replicates the 
single particle EBK spectrum obtained from torus quan- 
tisation more precisely than it duplicates the exact single 
particle quantum spectrum. After inspection of Figs. ^ 
and 1^, we notice the same effect in the two particle spec- 
trum. This property of the trace formulas also accounts 
for the unresolved peaks in the semiclassical spectrum. 
When the spacing of two levels of the EBK spectrum is 
very small, our truncated trace formulas may not resolve 
them, regardless of the spacing of the corresponding lev- 
els in the quantum spectrum (cf. FigJJ). 

Comparing Figs. [| and ||, we observe generally good 
agreement between the quantum and semiclassical spec- 
tra. Still, there are some apparent inconsistencies, for 
example, the two tall peaks in Fig. u. These are the 
two sets of unresolved levels in Fig. U in each case an 
octet and a quartet. The reason for the discrepancy is 
the level spacings are smaller than the smoothing width 
a, in contradiction to the assumption above, so that the 
peak height does not equal the degeneracy. In fact, the 
peak heights observed are rather close to 12 since the 
octets and quartets are very nearly degenerate on the 
scale of a and act almost like a 12-fold degenerate set. 



9 



It is not perfectly 12 due to the fact that the degeneracy 
is not perfect. However, we also observe that the inte- 
grated weight under the peak is consistent with a set of 
12 energy levels. Other inconsistencies in Fig. |^ occur for 
the same reason. Of course, overall improvements can be 
made by including more orbits. As well, the artifacts of 
the single particle spectrum (some examples are marked 
by an "X" in Fig. |) which arise from errors in the cross 
term, presumably decrease when corrections to the sin- 
gle particle trace formula are incorporated into the cross 
term. These preliminary numerical findings support our 
analytical results, which we now test in the rather differ- 
ent context of a chaotic billiard. 



VI. TWO PARTICLE CARDIOID BILLIARD 



Gutzwiller amplitudes are computed using the standard 
prescription for the stability of free flight billiards (see, 
for example, 0.) 



*2a: 5.195 3a: 6.584 4a: 7.102 '4b: 9.237 5a: 10.38 






:: 10.94 7a: 11.26 '6b: 11.84 5b: 11.98 '8b: 13.33 







'10b: 14.21 6c: 14.22 7b: 15.26 7c: 17.22 '8c: 19.09 




In this section, we study the problem of two identical 
noninteracting particles evolving in the cardioid billiard, 
which is fully chaotic |^^. Since the billiard has a re- 
flection symmetry, all the quantum states are either even 
or odd (this symmetry should not be confused with the 
symmetric/antisymmetric symmetry due to particle ex- 
change.) In the subsequent analysis, we will exclusively 
use the odd spectrum. The reason for this is to avoid the 
additional complication of diffractive orbits which strike 
the vertex. Classically, these orbits are undefined and are 
therefore not included in the standard Gutzwiller the- 
ory. Studies of diffractive effects in trace formulas can be 
found in [ ^8|j2g[| . The latter reference explores the specific 
application to the cardioid and shows that diffractive or- 
bits are important in describing the even spectrum but 
are largely absent from the odd spectrum. 

We could proceed as before by doing an explicit semi- 
classical analysis of each term in the decomposition of the 
two particle semiclassical density of states ( ll|) . However, 
we can simplify the analysis by removing single particle 
dynamics from the discussion. That is, we will focus ex- 
clusively on those quantum mechanical and semiclassical 
quantities that inherently describe two particle dynam- 
ics. More specifically, we compare the Fourier transform 
of the dynamical term 



Fr(L)=.F{/5i*pi(fc)} 



(43) 



with its quantum mechanical analogue which we define 
to be 



Fr{L)=T{p2{k) 



pi* pi{k) ~2pi* pi{k)}. (44) 



The integral operator J-' will be defined precisely below. 
In the semiclassical transform ( p3| ) , we use (|3^) expressed 
in terms of the wavenumber. Here, Fi and F2 are peri- 
odic orbits in the fundamental domain (i.e. the half- 
cardioid.) Orbit properties are discussed in [|30y2g| ] and 
some representative orbits are shown in Fig. 1^. The sta- 
bility matrices in the denominators of the single particle 



FIG. 6. Some of the shorter periodic orbits of the cardioid 
in the full domain. The label of each orbit includes the num- 
ber of reflections and also a letter index to further distinguish 
it. The two orbits *8b and *10b reflect specularly near the 
cusp, contrary to appearances while the orbit 4a misses the 
cusp. From |29| . 

In the quantum mechanical analogue (^), P2{k) is 
the quantum two particle density of states P2{k) — 
6 {k — kj) where the superindex / denotes the pair 

of integers (i, j) and kj — ^ kf + fc|. In (^), we sub- 
tract the smooth average part and the part which con- 
tains single particle dynamics. Using A = 37r/4, £ = 6, 
and K. = 3/16 in (|23|), the smooth term is 



/ 9 9 
V167r ^ 128 



' 256 167r 



The cross term is given by the general expressions (|2 
and (|2~ 



(45) 



A. Numerics for the Unsymmetrised Cardioid 

As before, we take a = 1 and use a standard sized car- 
dioid as in 1^ to obtain the single particle spectrum. 
In this section, we numerically compare the two particle 
quantum mechanics with the two particle semiclassics. 
We do this by making a direct comparison of the Fourier 
transforms in the reciprocal space of orbit lengths, L. In 
this space, we expect peaks at lengths which correspond 
to the Euclidean lengths of the full periodic orbits of the 
two particle system. For instance, if the full orbit F is 
comprised of particle 1 travelling on the orbit Fi and par- 
ticle 2 traversing a distinct orbit F2, we expect a peak 

Lp^. In the event that both particles 
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are on the same orbit F, we expect a peak at v^ir- In 
this way, any peak in the two particle spectrum can be 
attributed to the dynamics of a particular periodic orbit 
of the full classical phase space. 

We construct the two particle spectrum by adding the 
energies of the single particle spectrum. We include the 
first 1250 single particle energies which allows us to con- 
struct the first 766794 two particle energy levels repre- 
senting all two particle energies less than 6.8856 x 10'^. 




9 10 11 



FIG. 7. The Fourier transform of the dynamical (purely 
oscillating) part of the two particle density of states. The 
solid line is the transform of the quantum mechanical two 
particle spectrum and the dashed-dotted line is the 

transform of the semiclassical two particle trace formula 
(p3[). All relevant geometrical periodic orbits with length 
Lr = ^I/pj^ + Lp^ < 11 have been included. (Symbols de- 
scribed in the text.) 



In Fig. ^ we observe good agreement between the 
quantum and semiclassical results for L < 6.5 and L > 
10.3. In the region 6.5 < L < 10.3, there are appreciable 
discrepancies for the following reason. Recall that the 
amplitudes of the two particle trace formula (|3^) apply 
only to billiard systems whose single particle periodic or- 
bits are isolated. In the single particle cardioid problem, 
there exist orbits which are not well isolated in phase 
space, in fact two geometric orbits and a diffractive orbit 
are sometimes very close in phase space. For example, 
the two geometric orbits 4a and *106 together with the 
similar looking diffractive orbit 4o' (not shown) ||2^]. In 
this event, the stationary phase approximation underly- 
ing the Gutzwiller formalism fails as does the argument 
in that diffractive orbits do not affect the odd spectrum. 
As a result, whenever a two particle orbit in the full space 
is comprised of one or both particles on one of these prob- 
lematic single particle periodic orbits, the resulting two 
particle amplitude is inaccurate. (There is recent work 
on uniform approximations to account for such effects 
|^2| , unfortunately it seems not to apply to the cardioid 
which has the additional curious feature that the bound- 
ary curvature is infinite at the vertex.) 



A 


Ti 


r2 


La 


1 


i(*2a) 
i(*2a) 


4a 


7.562 


2 


i(*10b) 


7.565 



TABLE I. A few of the periodic orbits which conspire to 
give trouble around L=7.6. 



For a precise numerical comparison, we define the 
Fourier transform 



/oo 
-oo 



ikL 



f(k) 



(46) 



as a function of the conjugate variable L. Here, w(fc) is 
the three term Blackman-Harris window function fSlll 




ko < k < kf 
otherwise 



(47) 



with (ao,oi,a2) = (0.42323,-0.49755,0.07922). We 
choose ko and kf so that the window function goes 
smoothly to zero at the first and last eigenvalues of the 
two particle spectrum. Numerical integration of (ES) and 
( p^ ) using this definition of is displayed in Fig. ^. In 
the semiclassical transform, a total of 100 periodic orbits 
including multiple repetitions were used. 
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We now consider some specific examples. Consider first 
the peak structures "o" and "*" . In this region, the single 
particle trace formula is erroneous and these errors 
propagate through to the cross term and inevitably to 
the quantum mechanical transform. We have also com- 
puted the cross term using quantum mechanics, that is, 
using pi — pi — pi in (Q) and confirmed that these dis- 
crepancies do not arise (cf. Figs. ||and^). Thus, these 
discrepancies are due to errors in the semiclassical ap- 
proximations of the cross term. 



80- 




FIG. 8. Same as figure |7| except that the cross term in 
(^) is computed using single particle quantum mechanics. 
(Symbols described in the text.) 



accounted for by inspection Table II. 



40 




L 

FIG. 9. The cross term for the cardioid billiard calculated 
using quantum mechanics (solid) and periodic orbit theory 
(dashed-dotted) . 

As the Table II and Fig. |^ show, there are many in- 
stances where both of the single particle orbits consti- 
tuting the full orbit are poorly isolated. In view of this, 
both single particle Gutzwiller amplitudes are incorrect 
making the product even worse. This accounts for the 
gross inconsistencies in this region of the reciprocal space. 
The other discrepancies can be accounted for in a similar 
manner. 



For the rest of the discussion, A refers to a particular 
periodic orbit family in the full phase space with each two 
particle orbit in this family comprised of the single parti- 
cle orbits Fi and V2 and Xa are the lengths of the orbits in 
each family. Next, consider the peak structure at L « 7.5 
There are two families of orbits, Ai and A2 that are 
responsible for these peaks. The underlying structures 
of these orbits are shown in Table 1. Bearing in mind 
the two single particle orbits F2 are not well isolated (cf. 
Fig. H), the Gutzwiller amplitude of each r2 is incorrect. 
Consequently, the two particle Gutzwiller amplitude will 
also be incorrect, as Fig. |^ demonstrates. Let us look 
at the next peak structure. Clearly, the quantum peak 
heights are underestimated at L « 8.1 (X). We account 
for this by recognizing the two particle orbit structure 
involves the single particle orbit ^(*8b) which is an orbit 
which passes close to the vertex. More specifically, the 
orbit family A = 3 is composed of single particle orbits 
Ti = i(*4b) and F2 = 5(*8b) and the lengths of these 
two particle orbits are La = 8.109. As a final illustration, 
we consider the region 9.6 < L < 10.3. In this neighbour- 
hood, the semiclassics are particularly bad. This can be 



A Fi Fa La 

4 3a 4a 9.684 

5 3a i(*10b) 9.687 

6 4a i(*8b) 9.740 

7 5(*8b) i(*10b) 9.742 

8 4a 4a 10.044 

9 4a i(*10b) 10.046 

10 i(*10b) iflOb) 10.048 



TABLE II. A few of the periodic orbits which conspire to 
give trouble around L=10. 
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B. Symmetry Decomposition 

In this section, we explore the symmetry decomposition 
of the two particle problem and in particular the compar- 
ison of the symmetrized two particle quantum mechan- 
ics with the corresponding two particle semiclassics. We 
start by defining the smooth and oscillating symmetry 
reduced densities of states from ( pT] ) 



1 / 1 / fc 

Ps/A{k) = - ( (pi * pi){k) ± -^pi ( 



and 



{pi * pi){k) ± -^pi 



k 

71 



(48) 



(49) 



While the second term in (|49|) is a single particle den- 
sity, in a future paper |Q we will demonstrate that this 
term describes the physical situation in which two parti- 
cles are traversing the same periodic orbit, with the same 
energy and are exactly half a period out of phase. It de- 
scribes the effect of particle exchange on the spectrum 
and for this reason affects the symmetric and antisym- 
metric spaces differently and is based on the theory of 
Ref. [Q. Therefore, this second term also belongs to 
the two-particle dynamical term and we identify ( ^ ) as 
being a purely dynamical term. We want to compare it 
with the corresponding term in the symmetrized quan- 
tum densities of states. Hence, in analogy with the pre- 
vious subsection, we compare 



Jidyn 
S/A 



(50) 



and 



F^y^iL) =T{ps/A{k) - Ps,A{k) ~ pi * pi(fc)} . (51) 

where Ps/A{k) is the quantum bosonic (S) or fermionic 
(A) density of states. 



C. Numerics for the Symmetrized Cardioid 

In this section, we numerically compare the symmetrized 
quantum mechanics with the corresponding semiclassical 
quantities. In particular, we compute the transforms ( |50| ) 
and (pll) p4[|. The symmetrized quantum densities are 



(fc) - ^ 5 (fc - Jkf + kj\ +Y.5[k- V2kA , 



PA 



{k) = Y.s[k-^kf + k]) 



(52) 



using the same constraint on the energies as above. Of 
course, the sum of these symmetrized densities is the to- 
tal density of states. 



Before presenting our numerical results, we describe 
what we expect. First, all the peaks of the unsym- 
metrised two particle density should be present. In ad- 
dition, for each periodic orbit F, there should also be 
peaks at lengths Lp/V^ arising from the oscillating part 
of the single particle density of states. The results are 
shown in Fig. For the two particle density term of 
(pO|), we used the same 100 two particle orbits of section 
VIA while in the single density term we used all single 



particle orbits with length L < 11. As well, we included 
the single particle orbit i(*10h) (not shown in Fig. |^) 
which has a length L = 10.477. Fig.^ displays the peak 
structure in the reciprocal space up to L — 6.75. 




FIG. 10. The Fourier transform of the quantum and semi- 
classical symmetrized densities of states; the top is bosonic 
and the bottom is fermionic. In both cases the solid line 
is the transform of the quantum density of states and the 
dashed-dotted line is the transform of its semiclassical ap- 
proximation. Peaks with symbols are described in the text. 

We notice that most of the amplitude divides evenly 
between the symmetric and the antisymmetric densities. 
Nonetheless, there are exceptions such as the peak at 
L « 3.6 (o). Here, both terms of ( p9| ) contribute and the 
difference in the sign of the second term accounts for the 
uneven amplitude division between the two symmetrized 
densities. Semiclassically, we account for the peak struc- 
ture by noting that two different physical situations are 
responsible for the peak structure "o". First, there is 
the situation in which both particles are on the orbit 
F = ^(*2a) with no restrictions on the time phase differ- 
ence between the two particles. This contribution comes 
from the two particle density term resulting in a peak at 
a length v^Lr = 3.673 and produces identical structures 
in both densities. The second situation occurs when both 
particles are on the orbit F = (*2a) exactly half a period 
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out of phase. This contribution comes from the single 
density term at Ly-/^/2 = 3.673 and is explained more 
fully in 1^ . Since the second contribution comes with a 
different sign in the two symmetries, the amplitudes are 
different for the symmetric and antisymmetric spaces. In 
this particular case, it is stronger in the symmetric den- 
sity and weaker in the antisymmetric density, although 
the opposite may be true in other cases. 

In closing, we remark that the overall agreement be- 
tween the quantum and semiclassical calculations is good. 
The poorly reproduced peak just above L = 5 {+) comes 
from the single density term. This is just the poorly re- 
produced peak of the single particle density at L w 7 
shifted down by a factor of a/2. 



VII. CONCLUSION 

Initially, we developed a semiclassical formalism to de- 
scribe the two particle density of states. After deriving a 
trace formula describing two particle dynamics, we inves- 
tigated its structure and noted intuitive properties such 
as the additivity of the actions and topological phase fac- 
tors. As well, we briefly explained the structure of the 
full two particle orbits which come in degenerate families. 
As a first application, we wrote down a two particle trace 
formula for two identical particles in a billiard. The semi- 
classical symmetry decomposition involved formal sub- 
stitution of the semiclassical quantities into the quantum 
mechanical expressions for the symmetrized densities. In 
a future paper pO[ | , we show how these formal expressions 
emerge directly from the classical structures. 

Following these general considerations, we studied two 
identical noninteracting particles in a disk and in a car- 
dioid. In each case, we find that the formalism correctly 
reproduced the full and symmetrized densities of states. 
In the integrable problem, we found that our formalism 
replicates the two body EBK spectrum more precisely 
than the quantum spectrum, suggesting a deep connec- 
tion between periodic orbit theory and EBK quantization 
for integrable systems. In the chaotic cardioid billiard, 
we note that the single particle orbits which pass close 
to the vertex lead to inconsistencies in the Fourier trans- 
form of the semiclassical density of states. Clearly, our 
formalism fails here because the Gutzwiller theory itself 
fails for these "semi-diffractive" orbits. For all other or- 
bits, the two particle trace formula works very well. 

The techniques employed here involve the classical 
phase space of each particle. In a future paper pc|], we 
derive the same results by working in the full two par- 
ticle phase space. This approach has the advantage of 
being more general than what we have presented here. 
Nonetheless, it is conceptually useful to see how the same 
structure emerges from these two distinct points of view. 
We would also like to incorporate interactions between 
the particles. Such a project would undoubtedly require 



working in the full phase space since it is no longer true 
that the full density of states is the convolution of the sin- 
gle particle level densities. This provides an additional 
motivation for working out the noninteracting problem 
in the full phase space as a first step towards the more 
ambitious goal. This full phase space analysis also gener- 
alises more readily to more particles. Finally, it has the 
conceptual advantage that the spurious endpoint contri- 
butions discussed in [VD| and Appendix B do not arise 
and therefore need not be explained away. 

It may be argued that interacting many body sys- 
tems are too complex to be accessible to the semiclas- 
sical method. However, given the intractability of the 
many body problem, there may well be questions which 
semiclassical theory can answer. In particular, we have 
in mind the applications of semiclassical theory to meso- 
scopic physics [|l9|. Here, our seemingly academic study 
of billiard systems finds physical applications in the con- 
text of nanostructures. For example, the disk billiard 
can serve as a realistic lowest-order approximation to the 
mean field of the electrons in a circular quantum dot . 
In fact, many phenomena in ballistic mesoscopic systems 
can, at least qualitatively, be described by using quantum 
billiards with independent particles as physical models. 
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VIII. APPENDIX A: NONIDENTICAL 
PARTICLES 

As we have mentioned, most of the discussion still ap- 
plies if the two particles are not identical. Another sit- 
uation is a single particle in a separable potential. For 
example, in two dimensions, one could have V{x, y) = 
Va{x) + Vb{y) in which case, the dynamics in the x direc- 
tion are completely uncoupled from the dynamics in the 
y direction so that the system is formally the same as if 
there were distinct particles executing the x and y mo- 
tions. The formalism presented above follows in a natural 
way. The main differences are that one no longer consid- 
ers the symmetrised density of states since the symmetry 
of particle exchange no longer exists and secondly there 
are two distinct cross terms so that ([ll|) is replaced by 



P2[E) = pia * Plb{E) + pia * Plb{E) 
+ Pia * Plb{E) + pia * Plb{E), 



(53) 



where the indices a and b refer to the two distinct par- 
ticles, while the indices 1 and 2 still refer to one or two 
particle densities of states. 
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Imagine, for example, that we have two nonidentical 
particles in distinct billiard enclosures. We introduce 



two parameters, aa = 2ma/h and ab 
smooth term is replaced by 



2mb/h. The 



In the special case of identical particles, it is simple to 
check that this expression reduces to ( |30| ) . For lack of an 
immediate physical context, we do not explore this case 
any further. 



Pla * Pib{E) 



167r2 

.1/2 



E 



la 

1/2 1/2 

647r 



(gg + ab)AK. 
An 



(54) 



The cross terms each separately have the same structure 
as the cross term for identical particles. Obviously, they 
are no longer equal to each other, but functionally lit- 
tle has changed. It is just a question of inserting the 
relevant information from the different smooth and oscil- 
lating densities of states of the two particles. Following 
the same logic as before, we find 



IdE) 
Ik{E) 



> cos 



4^' jr^Trb 

1 b 



2^ -7^ COS ( 



87r3/2V^ 



Trb 



^ArfcCos($rJ, 



(55) 



where $ri, = VObELr^ - aTi,n/2. For pu * pia{E), we 
just interchange a and b. 

The formula for pia * Pib{E) still has the same ba- 
sic structure, but should obviously use the distinct peri- 
odic orbits for particles a and b. In particular, Eqs. ( p^ ) 
and (|l^) still apply, but with two important differences. 
Firstly, the double sums over periodic orbits are now la- 
belled by the distinct periodic orbits of the two parti- 
cles. Secondly, the energy partition will change due to 
differing masses. The criterion of stationary phase will 
still specify that the two particles have the same period, 
but relations such as (|2^) and (^8|) do not apply since 
they assume equal masses. The generalisations are rather 
straight-forward to determine. For example, the saddle 
energies (ESh are replaced by 



Ef) 



rubLy^ 



E - Eq 
E 



TUbLy^ 



IX. APPENDIX B: SPURIOUS END-POINT 
CONTRIBUTIONS FOR THE CARDIOID 

Here we evaluate the cross term integrals exactly for 
isolated periodic orbits. This allows us to do an asymp- 
totic expansion to explicitly demonstrate that the addi- 
tional endpoint contributions not included are spurious. 
We must evaluate the integral 



mgLp 



jTibL'i^ 



while the general dynamical expression for billiards ( |30| ) 
is replaced by 



Pla * Pib{E) 



(27r)3/2 



E 

COS 



imaEl 



rubLf 



-ArAEo)ArAE - Eo) 



n TT 
2 ~ 4 



(57) 



pi * pi{E) 



f dep,{e)p,{E ^ e) 
Jo 



(58) 



where pi(e) is given by the Weyl expansion (^ and pi{e) 
for a billiard with isolated orbits is given by 



^^(^)"2^E 

This gives 
pi*pi{E) w 



det 



(Mr - /) 



(59) 



det Mr - / 



a3/2^/, + ^1/2^/3) (60) 



87r2 



167r2 



2tt 



where 

h = 
h = 



de 







-^_cos(v/a(i^-e)Lr-ar- 



E 



1 1 



de- ^ , 

VeVE-e 

1 



cos 



{V^iE-e) 



- e)Lr - f^r-^ 



A — — ;= COS V aELr — (Tr — 
y/E y ' ' 2 



(61) 



If we evaluate the first two integrals exactly, we get 



(56) 



/i = 



a 



l/2Lr 



COS ( $r 



&r - - ) - cos 



(*r-^)) 



(62) 



and 



I2 = TT cos (j)r Jo ($r) — tt sin (j)r Hq ($r ) 
($r 



27r 



JuELy 
2 

'aELr 



■cos(0r) 



(63) 
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where $r = VaELr, (/r = — OT7i'/2, Jo is a zero-order 
Bessel function and Hp is a zero-order Struve function. 
In the second hne of (^3|), we have used the asymptotic 
expansions of these two functions. 

In both Ii and I2, we note that asymptotically there 
are terms with two distinct structures. The first are 
terms which are sinusoidal in ^/E and correspond ex- 
actly to what was used as the cross term for the cardioid 
{i.e. Eqs. ( |2^ ) and ^6|)). There are also terms which 
are nonsinusoidal in E. In Ji, this comes directly from 
the upper endpoint of the integral while in I2 it comes 
from the expansion of the Struve function. In each term, 
the nonsinusoidal terms arise from the endpoint around 
e = E which, as we argued in section [VP , corresponds 
to an unphysical situation. Therefore, keeping only the 
asymptotically appropriate term (i.e. the oscillatory one) 
yields the correct behaviour for the cross term. 

A similar analysis would yield similar results for the 
spurious endpoint contributions in the cross term of the 
disk billiard and the dynamical term of either billiard. 
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